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ABSTRACT 



Space missions to NEOs are being planned at all major space agencies, and recently 
a manned mission to an NEO was announced as a NASA goal. Efforts to find and select 
suitable targets (plus backup targets) are severely hampered by our lack of knowledge 
of the physical properties of dynamically favorable NEOs. In particular, current mission 
scenarios tend to favor primitive low-albedo objects. For the vast majority of NEOs the 
albedo is unknown. 

Here we report new constraints on the size and albedo of 65 NEOs with rendezvous 
deltaV < 7 km/s. Our results are based on thermal-IR f lux data obtained in the 



framework of our ongoing (2009-2011) ExploreNEOs survey ([Trilling et al.ll2010l ) using 
NASA's "Warm Spitzer" space telescope. As of 2010 July 14, we have results for 293 
objects in hand (including the 65 low-deltaV NEOs presented here); before the end 
of 2011 we expect to have measured the size and albedo of ~ 700 NEOs (including 
probably ~ 160 low-deltaV NEOs). 

While there are reasons to believe that primitive volatile-rich materials are univer- 
sally low in albedo, the converse need not be true: the orbital evolution of some dark 
objects likely has caused them to lose their volatiles by coming too close to the Sun. For 
all our targets, we give the clos est perihelion dista nce they are likely to have reached 



(using orbital integrations from lMarchi et al.ll2009l ) and corresponding upper limits on 



the past surface temperature. Low-deltaV objects for which both albedo and thermal 
history may suggest a primitive composition include (162998) 2001 SK162, (68372) 2001 
PM9, and (100085) 1992 UY4. 

Subject headings: infrared: planetary systems — minor planets, asteroids: general - 
radiation mechanisms: thermal — space vehicles — surveys 
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1. Introduction 

To date, two near-Earth objects (NEOs) have been targeted by space missions, both yielding 
a wealth of fascinating and groundbreaking insights into the past and current state of the Solar 
System: NASA's Near-Shoemaker mission went into orbit around its target (433) Eros in 2000 and 
landed on it in the following year; the Japanese mission Hayabusa arrived at (25143) Itokawa in 
2005 and scrutinized the NEO for a few months. In June 2010, Hayabusa succeeded in returning 
asteroid dust samples to Earth. 

Given the remarkable success of these missions it is perhaps not surprising that robotic NEO 
mission concepts are being considered at space agencies across the planet, including NASA's 
OSIRIS-REx, JAXA's H ayabusa 2, and ESA's Marco Polo (see, e.g jLauretta et alj|201~ol ; lokada et al 



2010 



Michel et al.ll2009l . for recent updates on these missions). In a speech i n April 2010, Pre sident 



Obama announced the goal of a manned space mission to an asteroid (see lAbell et al.l 12009 . for a 
corresponding NASA mission scenario including robotic precursor missions). 

Finding a suitable target asteroid is one of the challenging aspects of mission planning. Targets 
are tightly constrained in terms of their orbital dynamics and physical properties. Furthermore, 
launch windows are usually tight and the planning process long. Unforeseen delays due to techno- 
logical or financial problems risk eliminating the nominal target; it is therefore generally advisable 
to plan for contingency or backup targets. 



Dynamics, Av: As discussed by lShoemaker Helinl (119781 ). mission cost depends chiefly on the 
required amount of propellant, which follows from the total specific linear momentum, Av, that 
must be imparted on the spacecraft for it to reach the target orbit. Minimizing At; is therefore a 
top priority for practical reasons. It is worth emphasizing that there are a large number of NEOs 
that are reachable at a lower Av than that required to reach Mars. 

A realistic assessment of Av depends on the specific mission scenario and timing, and must be 
evaluated on a target-by-target basis. A customary first-order estimat e is the Av of a Hohman n 
transfer orbit, which is an analytic function of the orbital elements (see IShoemaker Helinl 1 19781 ). 
Thus, in a first target selection process objects with sufficiently small "Hohmann-Aw" are identified. 
Only those objects need to be studied in detail. For the purposes of this study, we will refer to the 
"Hohmann- Aw" as Av (without qualifiers). 



Physical properties: For the vast majority of NEOs, including low-Aw objects, nothing is known 
about their physical properties. Frequently, however, mission concepts require the target to be 
within a given size or mass range, e.g., in order to enable the spacecraft to orbit the target. 
Moreover, the science goals of some current mission scenarios require their target to be a "primitive" 
object, translating into constraints on their albedo and thermal history (see below). 
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Primitive objects: Some meteorites contain surprisingly pristine material that has suffered very 
little modification since the early days of the Solar System. Their asteroidal parent bodies are of 
particular interest for some NEO missions (especially since both NEAR-Shoemaker and Hayabusa 
targeted S-type asteroids, which have undergone significant processing). 

Judging from meteorite analogues, asteroids with very low albedo (geo metric albedo pv ^ 
7.5 %) are very likely to be "primitive" and vice-versa ([Fernandez et al.ll2005l ). A word of caution 
applies to NEOs, however: their relative proximity to the Sun can potentially cause their surfaces to 
heat up to the point that thermal surface alterations occur (de-volatilization, chemical reactions, 
etc.). There are hence two necessary conditions for an NEO to have a primitive surface: a low 
albedo and an orbital history that never brought the perihelion too close to the Sun. 



ExploreNEOs: This i s the third paper de s cribing results from the ongoing (2009-2011) Explore- 

The primary goal of this survey 



NEOs survey (following iTrilling et al.l 120101 : lHarris et al.l 120111 ) 
is to measure the size and albedo of ~ 700 NEOs based on observations with NASA's "Warm 



Spitzer" space telescope. ITrilling et al.l (120101) des c ribe t he goals and methods of the survey along 
with results for the first 101 NEOs; lHarris et al.l (120111 ) check the accuracy of the ExploreNEOs 
results against values published in the literature (where available) and find diameters to be typically 
consistent within 20%, albedos within 50%. 

As of 2010 July 14, we have data for 293 NEOs in hand including 65 objects with Av < 7 km/s. 
By the end of the survey, i.e. before the end of 2011, we expect to have measured ~ 160 low-Av 
objects. We chose to publish this first batch of results in order to alert the community to the 
existence of our growing database of characterized low-Af objects. 



Overview of this work: In Section [2] we describe our photometric data. Our modeling approach 
is described in Section [3l and resulting diameter and albedo estimates are presented in Section [H 
In Section [5] we study the thermal history of our targets. We discuss the implications of our results 
in Section [6] and summarize our conclusions in Section [7J 



2. Warm-Spitzer observations 



The observat ions reported herein use the post-cryog enic ( "warm" ) mod e of the IRAC camera 
(jFazio et al.ll2004l ) onboard the Spitzer Space Telescope (jWerner et al.l 12004 ) . Each NEO target is 
observed in the two photometric channels (channels 1 and 2) with central wavelengths of around 3.6 
and 4.5 //m, respectively. Observations are built up from frames that alternate repeatedly between 
the two channels, such that the resulting fluxes are quasi-simu l taneou s. Further details on our 
observation design and data reduction are given in ITrilling et al.l (|2010l ). 
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Table 1. Spitzer data sorted by Av 



Av 


Object 


Time 


H 




r 




A 




a 




f36 


f45 


(km/s) 




(U 1) 








(AU) 


(AU) 


(deg.) 


(mjy) 


(mjy) 


4.632 


(25143) Itokawa 


2010- 


■May-15 14:45:09 


19. 


.20 


1 


.018 





.052 


74. 


.56 


1. 


.767 ±0.041 


6.033 ±0.072 


4.755 


1996 XB27 


2010- 


•Jul-12 06:41:36 


21 


84 


1 


.121 





.192 


51 


.92 





.0178 ± 0.0040 


0.0428 ± 0.0061 


4.887 


(10302) 1989 ML 


2009- 


•Aug-21 03:09:35 


19 


.50 


1 


.100 





.152 


56. 


01 





.229 ±0.015 


0.560 ±0.022 


5.276 


(99799) 2002 LJ3 


2009- 


■Sep-19 22:38:44 


18. 


10 


1 


.238 


0. 


.354 


46. 


.18 





.125 ±0.015 


0.392 ±0.020 


5.280 


2001 CQ36 


2010- 


Apr-15 13:27:01 


22. 


45 


1 


.069 


0. 


.125 


55 


.50 





.0215 ± 0.0044 


0.0787 ± 0.0081 


5.302 


(52381) 1993 HA 


2009- 


-Nov-18 12:09:41 


20 


.20 


1 


.257 





.402 


46 


.59 





.0303 ± 0.0052 


0.150 ±0.011 


5.328 


2000 YF29 


2010- 


■Feb-20 04:43:34 


20 


16 


1 


.015 


0. 


.123 


83. 


03 





.131 ±0.012 


0.467 ±0.020 


5.391 


(1943) Antcros 


2009- 


■Sep-15 00:09:47 


15. 


75 


1 


.548 


0. 


.951 


40. 


.17 





.189 ±0.013 


0.582 ±0.022 


5.486 


(138911) 2001 AE2 


2009- 


Aug-13 14:53:36 


19. 


.10 


1 


.330 


0. 


483 


41. 


78 





.0244 ± 0.0047 


0.0888 ± 0.0086 


5.487 


2006 SY5 


2009- 


■Jul-30 18:18:27 


22 


08 


1 


.093 


0. 


.136 


54. 


07 





.0312 ± 0.0051 


0.1167 ±0.0097 


5.555 


(162416) 2000 EH26 


2010- 


Apr-29 07:39:49 


21. 


70 


1 


.125 


0. 


.217 


51. 


.21 





.0378 ± 0.0058 


0.1181 ± 0.0099 


5.565 


(162998) 2001 SK162 


2009- 


-Dec-17 23:35:37 


18. 


.00 


1 


.135 





.451 


63 


.68 





.815 ±0.028 


3.184 ±0.051 


5.653 


(68372) 2001 PM9 


2009- 


Aug-20 10:32:54 


18. 


90 


1 


.130 


0. 


204 


53. 


62 


3 


.928 ± 0.059 


19.61 ±0.14 


5.719 


(12923) Zephyr 


2010- 


•May-21 04:54:42 


16. 


10 


1 


.477 





.876 


41. 


.50 





.253 ±0.021 


0.453 ±0.023 


6.041 


(85938) 1999 DJ4 


2009- 


Aug-13 15:40:17 


18. 


60 


1 


.409 





.688 


43 


.10 





.0172 ± 0.0039 


0.0620 ± 0.0072 


6.069 


(433) Eros 


2009- 


■Aug-25 06:16:47 


11. 


16 


1 


.702 


1 


218 


36 


50 


13.99 ±0.11 


40.23 ±0.20 


6.070 


2000 GV147 


2009- 


•Oct-24 06:09:17 


19 


03 


1 


.064 





.336 


74 


.34 





.0998 ± 0.010 


0.337 ±0.017 


6.086 


2000 XK44 


2010- 


Jan-11 03:12:17 


17. 


73 


1 


.288 





.848 


51. 


.85 





.0285 ± 0.0051 


0.1040 ±0.0094 


6.106 


1993 RA 


2009- 


■Nov- 17 02:57:14 


18. 


91 


1 


.191 


0. 


.587 


59 


20 





.0165 ± 0.0039 


0.0544 ± 0.0068 


6.130 


(177614) 2004 HK33 


2009- 


-Jul-30 16:25:59 


17 


60 


1 


.057 


0. 


.095 


63 


88 


5 


.895 ±0.071 


22.92 ±0.14 


6.191 


1998 VO 


2009- 


-Dec-13 03:44:37 


20. 


37 


1 


.024 





.021 


76 


.62 


4 


.481 ±0.063 


15.31 ±0.11 


6.196 


2006 SV19 


2009- 


■Sep-15 16:08:48 


17. 


76 


1 


.166 


0. 


672 


60. 


97 





.109 ±0.010 


0.435 ±0.019 


6.236 


2005 JA22 


2009- 


■Nov-10 03:50:52 


18. 


.47 


1 


.296 


0. 


489 


46. 


66 





.0908 ± 0.0093 


0.319 ±0.017 


6.240 


(1627) Ivar 


2010- 


Jun-16 08:39:40 


13. 


20 


1 


.933 


1 


.249 


27 


.62 


1 


.264 ± 0.033 


2.670 ± 0.048 


6.279 


(87024) 2000 JS66 


2010- 


Jan-10 19:24:49 


18 


70 


1 


.104 





567 


65. 


.82 





.0163 ± 0.0023 


0.0427 ± 0.0061 


6.323 


(22099) 2000 EX106 


2010- 


■Jan-11 09:47:46 


18. 


00 


1 


.029 


0. 


246 


79. 


60 





.225 ±0.015 


0.842 ± 0.027 


6.364 


2003 SL5 


2009- 


■Sep-19 23:06:19 


19. 


14 


i 

i 


.114 


u. 


1 P.A 

. 104 


DO. 


oz 





.289 ±0.017 


1.052 ±0.029 


6.364 


(35107) 1991 VH 


2010- 


-Mar-20 17:04:49 


16 


90 


1 


.172 


0. 


.660 


58. 


.78 





.113 ±0.011 


0.477 ±0.021 


6.379 


(172974) 2005 YW55 


2010- 


■Mar-24 09:16:52 


19. 


30 


1 


.242 


0. 


386 


44. 


.41 





.0615 ± 0.0078 


0.171 ±0.013 


6.405 


(65679) 1989 UQ 


2009- 


■Oct-15 20:50:49 


19. 


40 


1 


.129 


0. 


.237 


58. 


.64 





.405 ± 0.061 


2.264 ±0.043 


6.431 


(159402) 1999 AP10 


2009- 


■Aug-31 22:54:55 


16. 


40 


1 


.292 


0. 


.838 


52. 


.32 





.090 ±0.010 


0.278 ±0.016 


6.491 


(143651) 2003 QO104 


2009- 


•Aug-21 02:40:17 


16. 


00 


1 


.402 


0. 


801 


45. 


93 





.276 ±0.042 


1.052 ±0.031 


6.507 


1989 AZ 


2010- 


■Feb-19 07:05:07 


19. 


49 


1 


.003 


0. 


.044 


93 


.49 


8 


.998 ± 0.090 


46.28 ±0.19 


6.512 


2006 W0127 


2009- 


-Nov-04 01:09:24 


16. 


18 


1 


.528 





.889 


40 


.23 





.102 ±0.011 


0.350 ±0.018 


6.516 


(140158) 2001 SX169 


2010- 


■Jan-26 04:30:24 


18 


30 


1 


.253 


0. 


430 


47. 


.17 





.092 ±0.010 


0.362 ±0.018 


6.526 


(100085) 1992 UY4 


2010- 


■Fcb-07 01:04:57 


17. 


80 


1 


.226 


0. 


832 


54. 


76 





.317 ±0.017 


1.694 ±0.037 


6.534 


(85839) 1998 Y04 


2010- 


Apr-05 10:55:52 


16. 


30 


1 


.326 





.783 


48. 


.87 





.179 ±0.014 


0.792 ± 0.026 


6.539 


(68359) 2001 OZ13 


2010- 


■May-26 02:33:35 


17. 


60 


1 


.520 


0. 


896 


39 


56 





.0243 ± 0.0079 


0.0443 ± 0.0078 


6.577 


(11398) 1998 YP11 


2010- 


■Apr-11 12:59:19 


16. 


30 


1 


.480 


0. 


700 


36 


03 





.240 ±0.016 


0.830 ±0.027 


6.608 


(155334) 2006 DZ169 


2009- 


-Dec-20 23:59:59 


17. 


10 


1 


.636 





.871 


32 


.65 





,0548 ± 0.0072 


0.148 ±0.011 


6.609 


(175706) 1996 FG3 


2010- 


-May-02 05:19:51 


18. 


20 


1 


.213 





.418 


51 


.02 





.853 ± 0.028 


4.391 ±0.062 


6.610 


(52760) 1998 ML14 


2010- 


■Jan-26 05:06:11 


17. 


50 


1 


.075 


0. 


484 


69 


33 





.127 ±0.012 


0.463 ± 0.020 


6.628 


(138947) 2001 BA40 


2010- 


Jan-10 17:16:01 


18. 


40 


1 


.173 


0. 


.721 


59 


.14 





.0135 ± 0.0035 


0.0592 ± 0.0071 


6.652 


2002 QE7 


2009- 


■Sep-23 04:42:21 


19. 


33 


1 


.217 


0. 


.329 


47 


.87 





.0662 ± 0.0080 


0.188 ±0.013 


6.652 


(5626) 1991 FE 


2009- 


■Aug-21 03:23:36 


14. 


70 


1 


.689 


0. 


972 


33 


.08 





.434 ± 0.020 


1.209 ±0.031 
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In this work, we restrict ourselves to objects observed on or before 2010 July 14 with a ren- 
dezvous Av < 7 km/s. Av values are taken from La nce Benner's online l i st of Av for all NEOs0 
which is calculated from the orbital elements and the IShoemaker &: Helinl (| 19781 ) formalism. 



In Table [TJ we present the measured in-band fluxes and the observing circumstances as taken 
from JPL's Horizons ephemeris server. H magnitudes are assumed to be uncertain by 0.5 mag 
(see Section [3 ] ) . Ob servations carried out on or before 2009 November 4 have been presented in 



Trilling et al.l (|2010l ). later observations are new here. 



3. Thermal modeling 

We use an updated version of the t hermal-modeling pipeline used in prev ious publications 
resulting from the ExploreNEOs survey ([Trilling et al.l 120101 ; lHarris et al.l 120111) . For complete- 
ness, we briefly summarize the more detailed description given in iTrilling et al.l (|2010l ). Section 
13.11 presents the updates relative to the previous pipeline, chiefly an estimation of the statistical 
uncertainty of our results. Due to said update, s ome of our results d iffer slightly (but within the 
error bars) from the preliminary results given in ITrilling et al.l (|2010j) . A reanalysis of our entire 
data set is deferred to a later work. 

NEO fluxes at Warm-Spitzer wavelengths have significant contributions from reflected sunlight 
and from thermal emission. We are interested in the latter in order to calculate the target size and 
albedo using a thermal model. Therefore, in a first step, we estima te the amount of reflected solar 
radiation us ing the method first described by Mueller et al. (|2007l ). then refined in Trilling et al 



(|2008l . I2OIO1 ) . Briefly, we calculate the expected V magnitude based on the observing circumstances 
given in Table Q] and extrapolate to 3.6 and 4.5 /im fluxes using published values of the solar flux at 
those wavelengths and the Sun's V magnitude. We also ass ume the spectral reflectivity at Warm 



Spitzer wavelength s to be 1.4 times that in the V band (see ITrilling et al 



2008 



Harris et al 



200 



9; 



Trilling et al.ll2010l ). In-band thermal flux equals total flux minus reflected flux. In rare cases, the 
calculated reflected flux exceeds the measured flux, leading to unphysical negative thermal fluxes in 
channel 1. In these cases (which we expect are due to inaccurate H magnitudes and/or lightcurve 
effects), we drop the channel-1 flux from the thermal analysis. 

In-band thermal fluxes are color corrected to take account of the spectral breadth of IRAC's 
filters and the significant difference between the spectral shape of asteroidal thermal emission and 
the stellar-like spectrum assumed in IRAC flux calibration; color-correction factors for the reflected 
solar component are negl igible. Color - corre ction factors are calculated for each target using the 
IRAC passbands gi yen bvlCarev et al.l (20101 ). the observing circumstances, and a suitable thermal 
model. As found by lMueller et al.l (|2007l ) . the dependence of color-correction factors on the physical 
properties of the asteroid can be neglected. 



: http : //echo . j pi .nasa. gov/ ~lance/delta_v/delta_v. rendezvous .html 
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Table 1 — Continued 



Av 


Object 


Time 


H 




r 




A 


a 




f36 


f45 


(km/s) 




(UT) 






(AU) 


(AU) 


(deg.) 


(mjy) 


(mjy) 


6.661 


(66251) 1999 GJ2 


2010-Fcb-09 16:33:48 


17. 


.00 


1 


.751 


1.038 


30. 


.67 


0. 


0294 ± 0.0052 





.0512 ± 0.0066 


6.692 


(85990) 1999 JV6 


2010-Jul-12 14:35:39 


20 


00 


1 


.062 


0.116 


62 


.78 


0. 


978 ± 0.029 


4 


.922 ± 0.064 


6.703 


1998 SE36 


2010-Apr-28 10:35:18 


19. 


32 


1 


.294 


0.457 


42. 


.22 


0. 


0280 ± 0.0050 





1089 ± 0.0096 


6.738 


(5645) 1990 SP 


2010-Jun-15 07:08:46 


17. 


00 


1 


.780 


1.095 


30. 


80 


0. 


0496 ± 0.0071 





.237 ±0.015 


6.747 


(3671) Dionysus 


2010-Jun-12 23:01:07 


16. 


30 


1 


.136 


0.518 


62. 


60 


0. 


180 ±0.013 





.390 ±0.019 


6.750 


2002 HF8 


2009-Aug-24 23:44:37 


18. 


.27 


1 


.262 


0.449 


48 


.69 


0. 


1093 ± 0.0099 





.455 ± 0.020 


6.751 


(164202) 2004 EW 


2009-Aug-06 05:32:19 


20. 


80 


1 


.049 


0.158 


75. 


.37 


0. 


0459 ± 0.0062 





.135 ±0.010 


6.757 


2002 UN3 


2010-Feb-12 22:31:31 


18. 


60 


1 


.325 


0.503 


41. 


.88 


0. 


0343 ± 0.0058 





.0555 ± 0.0070 


6.791 


(90373) 2003 SZ219 


2009-Dec-20 00:33:14 


18. 


80 


1 


.313 


0.466 


42. 


.42 


0. 


.0332 ± 0.0057 





.0686 ± 0.0078 


6.817 


(40329) 1999 ML 


2009-Scp-03 00:25:51 


17. 


70 


1 


.344 


0.501 


41. 


.48 


0. 


197 ±0.014 





.623 ±0.023 


6.828 


(6239) Minos 


2010-Mar-20 17:55:44 


17. 


90 


1 


.122 


0.340 


61 


.27 


0. 


114 ±0.012 





.28 ±0.17 


6.830 


2005 EJ 


2010-Mar-19 13:32:33 


19. 


87 


1 


.225 


0.414 


49. 


.10 


0. 


0215 ± 0.0046 





.0530 ± 0.0067 


6.840 


(5646) 1990 TR 


2009-Dec-01 18:48:50 


14. 


30 


1 


.857 


1.679 


33 


.09 


0. 


0393 ± 0.0063 





.0809 ± 0.0085 


6.865 


2003 W07 


2009-Oct-09 15:38:14 


18. 


91 


1 


.237 


0.418 


50 


83 


0. 


225 ± 0.022 





.517 ±0.024 


6.867 


1999 RH33 


2010-Apr-04 02:25:30 


19. 


13 


1 


.286 


0.463 


43 


.65 


0. 


0177 ± 0.0042 





.0371 ± 0.0057 


6.963 


(10115) 1992 SK 


2010-Apr-13 09:53:12 


17. 


00 


1 


.170 


0.307 


50. 


.53 


0. 


549 ±0.024 


1 


.896 ± 0.042 


6.963 


2003 BT47 


2010-Apr-21 05:07:07 


17. 


42 


1 


.289 


0.588 


48. 


67 


0. 


174 ±0.013 





.670 ± 0.024 


6.974 


(5620) Jasonwhcelcr 


2009-Aug-25 07:00:34 


17. 


00 


1 


.321 


0.620 


48. 


68 


0. 


258 ±0.016 


1 


.252 ± 0.033 


6.981 


(152563) 1992 BF 


2010-Apr-03 23:55:21 


19. 


80 


1 


.125 


0.236 


53 


.50 


0. 


282 ±0.017 


1 


.377 ± 0.034 


6.994 


(138971) 2001 CB21 


2009-Dec-06 22:09:34 


18. 


10 


1 


.033 


0.272 


79 


.87 


0. 


176 ± 0.014 





.591 ± 0.023 



Note. — Times given are roughly mid-observation as measured on Spitzer. H is the absolute optical magnitude (taken from 
Horizons and assumed to be uncertain by ±0.5 mag), r and A are heliocentric and Spitzcr-centric distance, resp., a is the solar 
phase angle. f36 and f45 are the flux at ~ 3.6 and 4.5 fim, resp. (channels 1 and 2). 
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Diameter and albedo are estimated from the final thermal fluxes using the NEATM ([Harris 

19981 ). The NEATM contains a dimensionless parameter r] that describes the effective surface 



temperature. iTrilling et al.l (|2008l ) found that data quality does not usually allow to fit rj to Warm- 
Spitzer data of NEOs, but that reasonable estimates of diameter and albedo can still be obtained 
by assuming an em pirical linear r elati onship between r\ and solar phase angle a. That relationship 
was established by Delbo' et al.l (120031) : here we use an updated relationship (based on a slightly 



larger data set) by lWolters et al.l (|2008l ): 

7] = (0.91 ± 0.17) + (0.013 ± 0.004)a (in deg.). 



(1) 



3.1. Monte-Carlo approach 

In order to provide a realistic estimate of the uncertainty in our diameter and albedo results, we 
use a Monte- Carlo approach in which various source s of uncertainty ar e considered: the measured 



flux uncertainty, the calibration uncertainty of 5 % (jCarey et al.ll2010l ). the uncertainty in H (see 
below), and the NEATM temperature parameter rj (which is assumed to vary by ±0.3 around its 
nominal value; see below). For each observation, we generate 1000 sets of random synthetic fluxes 
normally distributed about the measured value and with a standard deviation equal to the root- 
sum-square of the measured flux uncertainty and the 5 % calibration uncertainty. Analogously, 
Gaussian distributions of H and r/ values are used in the fit. The distribution of albedo results 
(diameters to a lesser extent) is strongly non-Gaussian, see Fig. [U We hence adopt the median of 
our Monte-Carlo results as the nominal value and asymmetric error bars to encompass the central 
68.2 % of the results. Additionally, we determine the percentage of albedo results falling into albedo 
bins (see Section E|). 



AH: In the analysis of purely thermal observations of asteroids, H is know n to have a negligible 
influence on best-fit diameters but impacts p\ directly ( Harris &; HarriJl997 ). In our case, however, 
H is also used in correcting for reflected solar flux and hence influences the calculated thermal 
flux contribution. We therefore decided to vary H within the Monte-Carlo fit. The correction for 
reflected sunlight becomes more critical as more reflected sunlight is contained within the measured 
flux; it is hence particularly important for high-albedo objects and for objects observed at large 
heliocentric distance. 

Prop agating AH into ther mal-flux uncertainties is an update of the thermal-modeling pipeline 
relative to ITrilling et al.1 (|2010l ). While this does not change the calculated nominal thermal fluxes, 
it does change their uncertainties and hence the relative weight with which they enter the \ 2 
minimization procedure, leading to somewhat different diameter and albedo results. For practically 
all targets, however, the corresponding diameter change is < 1 % and hence negligible. 

For most of our targets, published values for H, let alone AH, are unavailable. An observing 
campaign to measure H for a number of our targets is currently underway; results will be reported 
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Fig. 1. — Histogram of the distribution of diameter and albedo values resulting from the Monte- 
Carlo procedure for 2001 SK162. Both histograms are overplotted with the respective best-fit 
Gaussian functions. Note that the diameter distribution is rather consistent with being normally 
distributed, while the albedo distribution is clearly not. 



in a separate paper. For the time being, we fall back to the appro ximative H values given by the 
Hori zons ephemeris server which are Notorious* -mttftlfl ^ E^^I B lp a -W «. ,1 
20081 ). For now, we adopt AH = 0.5 mag throughout, see Section [4.11 for exceptions. 



At?: Changes in t he assumed r? value lead to significant changes in diameter and albedo, see e.g. 
Harris et al.l (120111 ). The quoted uncertainties in Equation (pQ) lead to a final A ?? ~ 0.3 for a ~ 50°, 



a typi cal value for our sample. This uncertainty estimate is corroborated by iRvan &: Woodward 
(J20ipJ), who found a typical ry value of 1.07 ± 0.27; we caution, however, that their sample is 
dominated by large main-belt asteroids, whose thermal properties are rather distinct from our 



2 It should not be forgotten that Horizons is not designed to calculate H magnitudes but ephemerides, at which 
it does an excellent job. 
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sample of small NEOs. For our Monte-Carlo analysis, we therefore adopt a Gaussian distribution 
of 77 values scattering about the nominal value of 0.91 + 0.013a with a standard deviation of 0.3. 
Unphysical 77 values below 0.5 are discarded. 



4. Sizes and albedos 

Our diameter and albedo results are given in Table [2] along with their statistical uncertainties 
estimated from the Monte-Carlo analysis described above. In order to illustrate the implications 
of our albedo results on surface mineralogy, we also determine the probabilities Pi— P4 with which 
the albedo falls within one of four albedo bins; these probabilities are estimated as the fraction of 
Monte-Carlo albedos falling within the respective bin. The albedo bins are designed to correspond 
to taxonomic types as closely as possible, particularly for the purpose of identifying primitive NEOs: 



pv < 7.5 %: As shown bv lFernandez et al.l (|2005l ). albedos in this range strongly indicate a 
primitive surface composition 

7.5 % < py < 15 %: While objects in this albedo range are still likely to be primitive, some 
of them may be more akin to (silicate rich) S types 

15 % < py < 30 % Objects in this albedo range are most likely S or Q types (or M types, 
but those are relatively rare among NEOs) 

30 % < py More exotic compositions, e.g. E types. 



Table 2. Results of our Spitzer observations and of the dynamical analysis described in Section [5] 



Ad Object D py 

(km/s) (km) 

4.632 (25143) Itokawa 0.319+.°°,;S QA1 -°'.1S 

4.755 1996 XB27 - 084 io'oi2 °- 48 -0.19 

4.887 (10302) 1989 ML 0.248+.° oil 0.47+° ™ 

5.276 (99799) 2002 LJ3 0.503^0085 - 43 t°Tl8 

5.280 2001 CQ36 °- 068 io 012 °- 41 io'l6 

5.302 (52381) 1993 HA 0.337+.°° 9 .g 0.140+.°; Q 7° 



5.328 2000 YP29 0.244I" "*i 0.251 



+ 0.041 n +0.154 
0.038 U - ZD -0.095 



5.391 (1943) Antcros 2A8+°H 0.145t°; ^ 

5.486 (138911) 2001 AE2 0.352 +.°'°>£| °- 34+ .°i'l6 

5.487 2006 SY5 O.090+.°;°5| 0.34+°'i4 
5.555 (162416) 2000 EH26 °- 141 too24 °- 181 io'o73 
5.565 (162998) 2001 SK162 1.94+°-^ °- 031 io'.015 
5.653 (68372) 2001 PM9 1 - 73 t°'.il °- 0180 -o'.0080 
5.719 (12923) Zephyr i- 86 !?'.!! °- 21 -o'.ll 
6.041 (85938) 1999 DJ4 °- 478 to 082 °- 28 to'l3 

6.069 (433) Eros 30.7+^° 3 ° 0.065l°;°^| 

6.070 2000 GV147 0.502+°- J°* 0.185+.°; g4 
6.086 (217807) 2000 XK44 0.73+o Jg 

6.106 1993 HA °- 358 to 052 °- 40 io'l6 

6.130 (177614) 2004 HK33 0.94+°;J 7 °- 189+ .o'.081 

6.191 (192559) 1998 VO 0.216+.°,;°§2 °- 30 to'.13 

6.196 (212546) 2006 SV19 1.06+.°'^ °- 129 i<xa5S 

6.236 2005 JA22 0.67+°J 7 0.164+°- 

6.240 (1627) Ivar 9 - 4 l?'.3 °- 094 lS'o51 

6.279 (87024) 2000 JS66 0.312+°°|;j 0.63+°|3 

6.323 (22099) 2000 EX106 0.621+° ™ 0.29+°,J° 

6.364 2003 SL5 0.337+.°° 6 i 3 O^S+^H 

6.364 (35107) 1991 VH 1.12+° ^ - 27 lo'.i° 

6.379 (172974) 2005 YW55 0.342+°;°™ 0.30+%-^* 

6.405 (65679) 1989 UQ 0.73+.°J| 0.060+.°;°^° 

6.431 (159402) 1999 AP10 1.20+.°;^ 0.35+.°;^ 

6.491 (143651) 2003 QO104 2.29+° J* °- 137 lo'o61 

6.507 1989 AZ 1M tolg °- 025 to'.oll 

6.512 (218863) 2006 W0127 1.70+.°;3° 0.21±°-} 7 

6.516 (140158) 2001 SX169 0.58+.°; J° 0.26+.°}° 

6.526 (100085) 1992 UY4 2.60±°™ 0.020+°;°ig 

6.534 (85839) 1998 Y04 l.SltoH °- 174 lS'.086 
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0. 


000 





227 





.443 





270 





.213 


0. 
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4.1. Reanalysis with updated H and G values for select objects 

In the analysis above, we assume G = 0.15 and use H magnitudes from the Horizons ephemeris 
service. While both assumptions are known to be problematic, we are constrained to use them in 
the 'mass production' of diameters and albedo for practical reasons: we are not aware of a central 
database of published H and G values; rather, values for each target have to be searched in the 
literature and are unavailable for a large majority. In order to minimize the induced uncertainties, 
we reanalyze our data for objects with published H and/or G values; see Table [3l We focus 
on objects which we find to have low albedo, as well as Eros and Itokawa. Where available, we 
also include published determinations of the size, albedo, or taxonomic type. The results of this 
reanalysis are given in Tabled! 



Low-albedo objects: For (175706) 1996 FG3, (65679) 1989 UQ, and (100085) 1992 UY4, our 
albedo results are low, as expected based on their known taxonomic classification. In the case of 
1996 FG3, there is also an excellent quantitative match between our results and the ground-based 
measure ments quoted in Table [3] Our diameter result for 1992 UY4 is formally in agreement with 
that by IVolquardsen et al.l (|2007l ). provided that their value is assigned a realistic uncertainty of 
~ 15 % to include systematic uncertainties (which are not discussed by the authors of that paper). 



'Reality checks:' While the general agreement between Expl oreNEOs results and other pub 



lished diameter and albedo results is discussed by Harris et al.1 (|201ll ). we consider it useful to 



check some of our results for low-At> objects against published values. By comparing Tables [3] and 
IU we find excellent agreement in the cases of (433) Eros, (25143) Itokawa, (10302) 1989 ML, (1943) 
Anteros, and (1627) Ivar. 



is hard to rec- 



(3671) Dionysus: Our nominal albedo result for Dionysus, py = 0.55 
oncile with its taxonim ic classification as Cf, type, for which a low albedo would be expected. 



-0.21 
-0.17' 



Harris Davied (|1999l ) report UKIRT mid-IR observations of Dionysus from which they derive 
py = 0.35 or py = 0.61, depending on thermal model. Due to "inad equate signal-to-n oise" , they 
were unable to constrain rj from their data, like in our case. However, lHarris &: Daviea also report 
ISO data, from which rj can be constrained (if marginall y so) to be ~ 3.1, in the u p per range of plau- 



sible rj values. This results in py = 0.16 (quoted after lHarris Lagerros 



2002, 



Harris &; Davies 



reject that result in favor of another thermal model, leading to much higher albedo). We have 
therefore repeated the analysis of our Dionysus data assuming rj = 3 (the phase angle of our obser- 
vations, a = 62.6°, is very similar to that of the ISO observations of lHarris Daviesl . 57.7°). The 
resulting albedo, py = 0.178jl ' 06 g, is in good agreement with Harris' NEATM result (py = 0.16) 
and consistent with a Cf, classification given the error bars. 



We note that bina ry NEOs including Dionysus were recently found to generally display higher- 
than-average rj values (jDelbo' et al.ll201ll ). probably due to regolith loss during binary formation. 



Table 2 — Continued 



Av Object D 

(km/s) (km) 



6.539 (68359) 2001 OZ13 0.62+°'^ 0.42+° 

6.577 (11398) 1998 YP11 1.73+_°;^ ' 176 too90 

6.608 (155334) 2006 DZ169 I.IsIqJI °-l 95 tu'.095 

6.609 (175706) 1996 FG3 1-9o1q'44 °' 026 ±0.012 

6.610 (52760) 1998 ML14 0.81+.°;^ 0.27+°'^ 
6.628 (138947) 2001 BA40 0.440+g;jJ|{; 0A2 -ois 
6.652 2002 QE7 0.32olJ5"°56 0.34l°' 14 
6.652 (5626) 1991 FE 3.96+J-g2 0.152+°J^ 
6.661 (66251) 1999 GJ2 0.90±°-^ 0.371°'™ 
6.692 (85990) 1999 JV6 °- 498 lo'o88 °' 076 -0.035 
6.703 1998 SE36 0.3431°°!^ 0.30l°, 
6.738 (5645) 1990 SP 2.20jl°;^ 0.062+° °™ 
6.747 (3671) Dionysus 0.89+.°,'" 0- 67 io 24 

6.750 2002 HF8 0.711°'}° 0.181+2 

6.751 (164202) 2004 EW 0.157+°'°^ 0.36+°'*j 



:.757 2002 UN3 0.310""'' „.% 0.67 



+ 0.052 n K _+0.35 
0.047 U,D -0.26 
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6.791 (90373) 2003 SZ219 0.3061°;°*° 0.58l°'2° 

6.817 (40329) 1999 ML 0.96+2'J| 0.154+° ,!,^ 0.117 0.370 0.347 0.167 0.608 0.994 1.000 521 408 406 £^ 

6.828 (6239) Minos °- 474 lo'o91 °- 56 io'.27 °' 000 °' 040 °' 123 °' 837 °' 301 °- 487 °' 694 708 556 466 I 

6.830 2005 EJ 0.230l°'°3g 0.43^°^* 

6.840 (5646) 1990 TR 2.03+.°'^ 0.65+.°' *,f 

6.865 2003 W07 0.681°'*° 0.109j;° ° 7 j 

6.867 1999 RH33 0.228l°'°*2 0.73^° ^ 

6.963 2003 BT47 1.151°;** 0.147+.°';!,*° 

6.963 (10115) 1992 SK 0.90+°'™ 0.34^°'^ 

6.974 (5620) Jasonwhcclcr 1.771°'*° 0.0941°'°°° 

6.981 (152563) 1992 BF 0.5ll° " 0.084l°° 77 

6.994 (138971) 2001 CB21 0.5781°'*°° 0.24l°;*jj 



Note. — Pi— PA denote the probability of the albedo falling within each of the four albedo bins described in the text (bin boundaries are 0.075, 
0.15, and 0.3; primitive objects should display a high pi). The last six columns, q xx and T xx describe the orbital and thermal history (see 
Scction^J. E.g., 9io% denotes the minimum perihelion distance that the object reached to within a probability of 10%, T^Qe/ is the corresponding 
temperature to which the surface was heated. Note that some of the diameter and albedo results herein are superseded in Table l4l 
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Table 3. Published physical properties for objects considered in Sect. 14.11 



Object 


H 


G 


D(km) 


Pv 


Bin a 


Taxo b 


Reference 


(433) Eros c 


10.46 ± 0.10 


0.18 


23.6 a 


0.22 b 




S 


1,2 


(25143) Itokawa 


l».Oi_g 0g 


n 29+ 07 


0.327 ±0.006 


0.26 ±0.02 




S(IV) 


3,4 


(10302) 1989 ML 






0.28 ±0.05 


0.37 ±0.15 




E 


5 


(175706) 1996 FG3 


17.76 ± 0.03 


-0.07 


~ 1.9 


~ 0.04 


Y 


c 


6,7,8 


(1943) Anteros 


15.9 ±0.2 


0.23 


~ 2.2 


~ 0.16 




L 


9,10,11 


(65679) 1989 UQ 


19.5 ±0.3 










B 


11,12 


(162998) 2001 SK162 












T 


11 


(100085) 1992 UY4 


17.71 ±0.10 




1.7 


0.05 




P 


13,14 


(152563) 1992 BF 












x c 


15 


(12923) Zephyr 












S 


16 


(1627) Ivar 


12.87 ±0.10 




9.1 ± 1.4 


0.15 ±0.05 




s 


16,17 


(85990) 1999 JV6 














15 


(3671) Dionysus 


16.66 ± 0.30 




1.1-1.5 


0.16-0.31 


Y 


c„ 


15,18,19 



a Y indicates that the object is known to be binary 
b Taxonomic type 

sec iTrilling et all (hold) for a disc ussion of the values adopted for Eros, which was observed at a nearly pole-on 



viewing geometry, and lMuellerl l|2007l ) for a discussion of its G value 



Li et al 



( 20041); ^bernardi et all (120081); (4) Volu me-equivalent 



Demura et al 



(200 



i); (5) 



Mueller et al 



I2QQ3); M 



References. — (l) lHarris fc Daviesl (|l999l ); (2)_ 
diame ter calculated from the volume given by 

d2006h ; (7) Mueller et al. , in preparat i on; (8 ) iBinzel et al.1 (|200lD ; (9) Adopted after 1 5.82 ± 0.14(IWisniewski et al 



Pravec et al 



19971 ) and 15.96 ± 0.14 dPrayec et al.lll998l). bo t h give G = 0.23; (10) Quoted after lHarris et al.l <l201ll) original 



data from lVeeder et al.1 J 19891); (11) ginzel et al.1 (|2004al ); ( 12) bravec et all (|l998h : (13) I Warner et al.1 J2OO6I )— no H 



uncertainty is g iven by Warner et al., 0.1 mag seems appropriate (or slightly conservative) given the low scatter of 
their data; ( 14) IVolquardsen et al.l (|2007l ) — -note that the i r erro r bars (not quoted herein) reflec t only the sta t istica l 



uncertainties, but not the systematics; fl5) lBus fc Binzell (|2002l ); (16) lBinzel et al.l (|2004bl ); (17) iDelbo' et all (|2003h 
determined D and pv of Ivar based on Keck mid-IR photometry, we here assume a 15% uncer tainty in D and 30% 
in pv as is usual for radiometric diameters. H is from Pravec et al. (unp ublished), quoted after 



(18) D ifferent, model-dependent, diameter and albedo values are given by Harris fc Daviesl (|1999l ); lHarris fc Lagerros 



Delbo' et al 



1 20031) 



(2002) — we quote the range of their adopted results. Harris fc Daviesl 1 19991) quote Pravec (unpublis hed) for the H 
given herein, no uncertainty value is stated, we assume a conservative uncertainty of 0.3 mag. (19) iMottola et al 
(|l997h 
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This may be expected to reduce the accuracy of our results for binaries in general. However, in 
the case of the only other known low-albedo binary, 1996 FG3 (see above), our results are quite 
consistent with those obtained otherwise. 



5. Thermal history 



It is well known that implantation of solar wind ions (Hapke 200 l|) and bomb ardment by 



micrometeorites can alter the spectroscopic properties of asteroids (jSasaki et al.ll200ll ). However, 



these aging processes affect only the topmost microns of the surface. This is not a problem for 
a sample collection experiment: current-technology sampling devices can sample material from a 
depth of a few centimeters, thus excavating below the space-weathered surface. 



However, iMarchi et al.l (|2009l ) have shown that the surfaces of a significant number of NEOs 
were heated by the Sun to very high temperatures that could induce surface alterations on previously 
primitive objects. Due to thermal conduction, a thin layer beneath the surface can be heat ed to 
similarly large temperatures; for typical thermal properties (jMuellerl 120071 : iDelbo' et al.1 120071) and 



rotat ion periods, the penetration depth of the heat wave is of the order of centimeters (ISpencer et al 
1989j), comparable to the digging depth of sample-taking devices. 



There is no clear-cut way to determine the maximum temperature to which an asteroid can 
be heated and remain in primitive condition. From laboratory studies of carbonaceous chondrite 
meteorites, to which primitive asteroids are generally believed to be related, we know different alter- 
ation processes that are characterized by different threshold temperatures (e.g. thermal breakup of 
organic macromolecules) . For instance, it has been determined by laboratory heating experiments 
that at 370 K the insoluble organic matter of the carbonaceou s meteorite Murchinson is degraded 
(the aliphatic C-H bond is lost) in approximately 200 years (jKebukawa et al.ll2010l ). The same 
authors also showed that the bulk organics of Murchinson are lost in only one year at 370 K (or 
200 years at 300 K). It is also known that the macromolecular phase in carbonaceous meteorites 
has a structure similar to refractory kerogen. The latter starts to break up - with production of 
oil and gas - when he ated above 420 K (I. Franchi 2008, personal communication). Furthermore, 
Lauretta et al.1 (|200ll ) have shown experimentally that volatile components (such as Hg) are re- 
leased from some CM and CV carbonaceous chondrites when the latter are heated above ~ 470 
K. 

The maximum temperature attained by an NEO is a function of the perihelion distance q. 
NEO orbits can evolve rapidly (see, e.g., iMichel fc Froeschlel Il997l ). hence the present q is not 
particularly indicative of the minimum q attained within the chaotic dynamical history. Hence, if 
it is a goal to send a spacecraft to a primitive object, the dynamical and thermal history of the 
target must be taken into account. 

An upper limit on the subsurface temperature is the surface temperature at local noon. Because 
spin axes of NEOs evolve on relatively short timescales (e.g. due to YORP, planetary encounters, 
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and possibly due to spin-orbit coupling) the whole surface is likely to have been subjected to 
temperatures (nearly) as high as that of the subsolar point, T$s- We calculate T$s as a function 
of q using the NEATM (see Section [3]), assuming the nominal albedo resulting from our Spitzer 
observations and 77 = 0.91. That latter value follows from Equation ([1]) for a = 0, hence providing 
an upper limit on temperature. 

With this in mind, the maximum surface temperature attained follows from the minimum 
q reached in the past. Due to the chaotic nature of NEO orbi ts, this questi o n mu st be treated 
probabilistica lly. We here use the orb ital evolution model by iMarchi et all (|2009i ). which was 
derived from iBottke et al.l (|2000i . |2002| ). For each of our targets, we extract from the work by 
Marchi et al.l the probabilities that the object reached a perihelion distance smaller than a grid 
of q values. Through interpolation, we determine q^o%, i-e. the q value which was reached with a 
probability of 50%. We call the corresponding subsolar temperature T 50 %. The odds are 50% that 
an object ventured to within q^o% of the Sun and that its surface was heated to T 50 % or more. We 
repeat this exercise for probability values of 10 and 90%; results are given in Table [2j 

Note, for instance, that (65679) 1989 UQ certainly appe ars primitive judgin g from its low 
albedo (this work) and its spectroscopic classification as B type (jBinzel et al.ll2004al ). 1989 UQ was 
considered repeatedly as a target of sample-return mission concepts to a primitive asteroid. How- 
ever, we show that this object has likely been heated to the point that its organic macromolecular 
matter has been broken up (T 50 y = 564 K — cf. Table [2]). The same applies to the binary system 
(175706) 1996 FG3. 

However, the final determination whether an object is to be considered primitive or not is 
beyond our scope. Rather, our aim is to provide the required input data for that determination. 
Given that objects may have exceeded some but not all of the threshold temperatures discussed 
above, the exact definition of 'primitive' depends on the scientific purpose at hand. Also, it is 
incumbent on the mission-design team to quantify the risk they are willing to take (note the 
probabilistic nature of our temperature determinations due to the chaotic orbital history of NEOs). 



Discussion 



6.1. Updated ExploreNEOs thermal-modeling pipeline 



For a fraction of our t argets, our data hay e been published pr eviously along with a straight- 
forward NEATM analysis ([Trilling et al.1 120101 ; lHarris et all l201ll ) . Due to the updates in the 
thermal-modeling pipeline presented in Section \3.1\ our results given in Table [2] supersede previous 
values where available. This difference, however, is always comfortably within the quoted error 
bars and too small to matter practically. A new analysis of the entire data set including new 
observations will be presented in a later paper. 



The new pipeline provides estimates of the statistical uncertainty of diameter and albedo 
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results. The uncertainties are distributed in a highly non-Gaussian way, especially for py, hence 
asymmetric "lcr" error bars are given. Additionally, we provide probabilities of py to fall into 
specific bins, which allows for more straightforward constraints on the taxonomic type and hence 
surface mineralogy. 



6.2. Potential spacecraft targets 

1996 XB27 and 1989 ML: These bodies are, along with the Hayabusa target Itokawa, the only 
objects with Av < 5 km/s within our sample. We find both objects to be very high in albedo , 
indicative of taxonomic types such as E (1989 ML was found to be E type bv lMueller et al.l 120071 ). 
The chances of them being primitive are very small. With a diameter of only 84^]J m > 1996 XB27 
is among the smallest celestial objects with known size. 



Primitive objects: Some of our albedo results suggest a primitive composition. We adopt a 
threshold value of 50% for p%, the probability of py < 7.5 %. A list of all measured physical 
properties of these objects is compiled in Table [5j There are significant object-to-object differences 
in thermal history. As discussed in Section [5j there is no clearly defined threshold temperature 
above which primitive material is metamorphosed. The least heated objects (at the 50% probability 
level) are 1992 UY4, 2001 SK162, and 2001 PM9. 1996 FG3 is the only known binary in the low-p v 
sample. 



7. Conclusions 

Of the 293 NEOs observed within the framework of our ongoing ExploreNEOs survey as of 
2010 July 14, 65 have Av < 7 km/s. Diameter and albedo measurements for the latter are presented 
in this work. Assuming that the rate of observations of low-Av objects stays as it is, the number 
of observed low-Av objects will increase to ~ 160 by the end of ExploreNEOs, i.e. before the end 
of 2011. Teams requiring a physical characterization of potential spacecraft targets are encouraged 
to contact us. 

Out of our 65 low-At> targets, 7 have low albedos indicating a primitive surface composition. 
These objects include a binary (1996 FG3) and three objects which stayed remarkably cool during 
their dynamical history, possibly cool enough to remain primitive: 1992 UY4, 2001 SK162, and 
2001 PM9. 



Michael Mueller gratefully acknowledges the Henri Poincare Fellowship, which is funded by the 
CNRS-INSU and the Conseil General des Alpes-Maritimes. The work of MM and MD is supported 
by ESA grant SSA-NEO-ESA-MEM-017/1. We thank Patrick Michel for helpful discussions. This 
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Table 4. Reanalysis of thermal data with updated H and G values (see Table [3]) 



Object 


D(km) 


PV 


Pi 


p2 


P3 


Pi 


D* 


pv* 


(433) Eros 


23.01 8 ; 3 


n 910+0.173 

u - zlo -0.099 


0.030 


0.250 


0.423 


0.297 


30.7+8°3° 


n65+ 075 

U.UDO_o 034 


(25143) Itokawa 


n 01 q+0.054 
U.OlO_ 44 


n 900+0.116 

U - ZOO -0.075 


0.000 


0.023 


0.550 


0.427 


0.319tg;»« 


41 +°- 20 


(10302) 1989 ML 


24n+ 043 


0.50l°;?| 


0.003 


0.010 


0.127 


0.860 


0.2481°;™ 


u -^' -0.19 


(175706) 1996 FG3 


1 S/I+ 056 
1.O4-0.47 


o.o42±g : g?f 


0.837 


0.140 


0.023 


0.000 


1 qn+ - 52 

l.yu_ 44 


0.026l°;° 29 


(1943) Anteros 


QC+0.72 
z -°°-0.59 


0.138±°:J°I 


0.150 


0.420 


0.367 


0.063 


9 /IO+0.69 


0-1451°;^ 


(65679) 1989 UQ 


72+ - 18 

u - ' z -0.14 


n nc -o+0.036 


0.753 


0.223 


0.023 


0.000 


n 7 o+0.18 
u - '•J-0.15 


nfin+ 059 

U.UOU_ .028 


(100085) 1992 UY4 


Z - OU -0.58 


0230+ 0190 

u.uzou_ .0090 


0.977 


0.020 


0.003 


0.000 


2 60+ - 70 


020+ 022 
u.uzu_ 010 


(1627) Ivar 


g 9+ 2 - 8 


0.128±°;J*i 


0.157 


0.447 


0.273 


0.123 


9.4+li 


n94 +0138 
U.U»1_o 051 


(3671) Dionysus 


oise+fi? 


u - oo -0.17 


0.000 


0.000 


0.047 


0.953 


0.891°;!! 


R7+ - 37 

u -°' -0.24 


Dionysus with 77 = 3 




n 17 q+o.092 

U.± 1 »_o.065 


0.020 


0.320 


0.553 


0.107 







Note. — In the last two columns, the diameter and albedo results for the 'nominal' H and G values are given for 
comparison (see Table [5J. In this reanalysis, 77 = 1.07 is assumed for (433) Eros (see Trilling et al.l [201ol ) . Two 77 



assumptions are made for (3671) Dionysus, the one used in the rest of this manuscript (see Equation |T}; upper line), 
and rj = 3 (lower line). For all objects except Eros and Dionysus, D and D* values are practically indistinguishable, 
implying that diameter results are not significantly impacted by the H uncertainty. 



Table 5. Potentially primitive objects as indicated by their low albedo (pi > 50 %) 



At; Object 


Pi 


D 


PV 




r 50% 


T 90% Taxo 


(km/s) 




(km) 




(K) 


(K) 


(K) 



5.565 


(162998) 2001 SK162 





92 


1 


9 ±0.4 





nq+0.03 
u,3 -0.02 


834 


462 


418 


T 


5.653 


(68372) 2001 PM9 





97 


1 


7 ±0.5 





m O+0.017 
u±o -0.008 


782 


488 


441 




6.405 


(65679) 1989 UQ 





75 





7 ±0.2 





nc;+0.04 
uo -0.02 


1082 


564 


410 


B 


6.507 


1989 AZ 





97 


1 


1 ±0.2 





n-3+0.02 
uo -0.01 


861 


560 


463 




6.526 


(100085) 1992 UY4 





98 


2 


°-0.6 





n9 o+0.019 
uzo -0.009 


498 


444 


417 


P 


6.609 


(175706) 1996 FG3 





84 


1 


0+O.6 
8 -0.5 





04 +0.04 
u ^-0.02 


743 


565 


482 


C 


6.738 


(5645) 1990 SP 





60 


2 


n+0.8 
z -0.7 





n fi +0.08 
UD -0.04 


743 


540 


456 





Note. — See Table [2] for the definition of p\. Where available, diameter and albedo results are 
from the reanalysis in Table HI otherwise from Table [2j As discussed in Section [5j the chances are 
10/50/90 % that the surface temperature has reached Tio/50/90% or above. 
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work is based in part on observations made with the Spitzer Space Telescope, which is operated by 
JPL/Caltech under a contract with NASA. Support for this work was provided by NASA through 
an award issued by JPL/Caltech. 
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